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ABSTRACT 

In this thesis the generation of equivalent networks by means of a 
congruent transformation using a variable transformation matrix applied 
to the admittance matrix of an n+1 node, common datum, network is 
studied. The ranges of the values of the elements of the transfor- 
mation matrix which yield passive RLC components for the equivalent 
network without mutual inductances are considered. When the trans- 
formation matrix is diagonal the solution involves a set of linear 
algebraic inequalities. In the general case, when all the elements 
of the transformation matrix are present, the solution involves non- 
linear algebraic inequalities. The region of passive component 
realizability consists of a hyper-volume in the -space of the elements 
of the transformation matrix, and the use of exponential searching 
functions to generate equivalent networks with design constraints 
is investigated. 
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Equivalent Networks - Historical Background 
1 . 1 Introduction 

It is well known that there is no unique solution to a synthesis prob- 
lem for realizing a given driving point or transfer function. A group of 
networks can always be constructed which have the same input immittance, 
or transfer function, or both. These groups of networks are called 
equivalent networks for the particular network function of interest. In 
general, equivalent networks can be generated by two distinct approaches . 

(a) Using different synthesis procedures for realizing the network 
function, such as Brune, Darlington, Bott-Duffin, Miyata or other pro- 
cedures for the case of input impedance synthesis. These procedures 
are well covered in the literature for input as well as transfer function 
synthesis, Ginllemin [3] , Van Valkenburg [8] . 

(b) Using a proper linear transformation of the admittance or im- 
pedance matrix of a given network such that the immittance between any 
desired node pair, or the transfer function relating the response in any 
node pair to the excitation at another pair, is kept invariant. 

The second method is a more general approach and is readily applic- 
able to computer design techniques . It leads to a technique for gener- 
ating a group of equivalent networks among which the best suited to the 
problem at hand can be chosen. 

The method of generating equivalent networks by linear transfor- 
mation was first introduced by Cauer [l] in 1929 in Germany and the 
first English paper on the subject was printed in 1931 by Howitt [2] . 
Guillemin [3] , discussing this subject, has shown that the same result 
as obtained by a transformation can be achieved by elementary oper- 
ations on the rows and columns of the network admittance or impedance 
matrix . Schoeffler [4] has shown that the transformation can also be 
applied to the branch matrices, and has also developed a technique of 
generating continuously equivalent networks [5] . 

In this chapter a discussion of the Cauer' s, Guillemin and 
Schoeffler' s method is presented. 
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1 .2 Equivalent Driving Point and Transfer Impedance 

The equilibrium conditions of a network containing n + 1 nodes 
with n ports defined by the terminals between each of the n nodes and 
datum, can be written in the following matrix form: 

[I ] = [Y ][V ] (1-1) 

where: 



[ V ] is the column matrix of the node to datum voltages 
[ I ] is the column matrix of source currents connected between each 
node and datum 

[ Y ] is the n x n node to datum admittance matrix of network 
The admittance matrix is formed in the following manner assuming 
no mutual coupling in the network. 

(1) Each diagonal element y is the sum of all the admittances 



connected to node i . 

(2) Each off diagonal element y i ^ j is the negative of the 
admittance connected between nodes i and j . 

The driving point impedance of a network excited at the input node 
i only is given by 



Z. 

l 



A.. 

li 



A 



(1-2) 



where: 

A = the determinent of the admittance matrix [ Y ] . 

th 

A. . = the minor obtained from Aby deleting the i row and 
the i^ column. 

Equation (1-2) can be obtained by expanding (1-1) and using Crammer's 
rule . 



It has been shown by Howitt [ 2 ] that if the admittance matrix 
[ Y ] undergoes a real non-singular transformation V 



L Y’] = [T] [ Y ] [r] T (1-3) 

where [F] is a square real non-singular matrix of the same order as 
[Y ] and T denotes transpose, then [Y‘] is the admittance matrix of a 
new network with the same number of nodes as the original network. 
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In general the input impedance at any node i and the transfer im- 
pedance between any two ports i and j of the new network is different from 

that of the original network. However Howitt [ 2 ] and Schoeffler [ 4 ] 

th 

have shown that if the i column of [ T ] consists of zeros except for 
th 

the i term which is unity^then the input impedance at port i remains 

★ 

invariant under transformation. If two such columns are chosen, the 
driving point impedance at the ports corresponding to the columns, as 
well as the transfer impedance between them, is invariant to the trans- 
formation. Example 1-1 illustrates the invariance of the input impedance 
under transformation (1-3) for a simple 3 node circuit. 

EXAMPLE 1-1 

Consider the resistance network of Fig. 1-la. 

1 




[Y ]= [G ] = 



Z. = 
in 



A. 

li 



2 

-1 

2 

3 



Using the following 

m= 1 8 



-1 

2 



transformation matrix: 

where a and b are arbitrary numbers yields 



[G 1 ] = 



2(l-a-2a 2 ) 

-b+2ab 



-b+2ab 




* 

An original proof based on matrix partioning is given in appendix I. 
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resulting in Fig. 1-lb where 

b(l-2a) 




A.. 

11 



in 



2b 

2b^(l-a+a^) - (-b+2ab)^ 



(a) Substituting a = 3, b = - will lead to the circuit of 

4/55 



in 




8/11 



(b) Substituting a = 3, b = -1 will lead to the circuit of Fig. 1-ld. 

1/15 



in 




1/3 



Fig. 1 . Id 



As can be seen the 3 networks are equivalent, in that their input im- 
pedance is the same. 

As can be seen from example (1-1), the choice of the elements of 
[T ] , other than the first column are arbitrary. If any of the elements 
in the transformed network are zero, the number of elements in the new 
network is different from the original network. 

If the network under consideration contains inductance and capaci- 
tances as well as resistances then we can write 

[Y ] = [G ] + -| [L _1 ] + S[C ] (1-4) 

where 

[ G ] = the network conductance matrix 

[ C ] = the network capacitance matrix 

[ L ■*■ ] = the network inverse inductance matrix 

By the distribution law of matrix multiplication [ 7 ] the transformation 
[ T ] can be applied to each of the above matrices to obtain the ad- 
mittance matrix of the new network . 

[ Y' ] = [ G' ] + ~ [ L ,_1 ] + S[C'] (1-5) 

Comparing (1-5) to (1-4) after tk§Jtt&ns.formation (1-3) 



[g 1 ] = [r] [g ] tr] T 


(l-6a) 


[c ]= [r][c ]Cr] T 


(1 — 6b) 


CL , ' 1 ] = [r][L' 1 ][r] T 


(l-6c) 



EXAMPLE 1-2 1 
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Consider the R, L, C network of Fig. l-2a, with the input impedance 



Z(s) = 



s 3 + S 2 + S 



3S 3 + 4S 2 + 3S + 1 



[G ] = 



3 

-1 



-1 

1 



[ L _ 1 ]= 



2 

-1 



-1 

1 



[C ] = 



0 

0 



0 

1 



Take [ T ] 



1 2 
0 -2 



; where 2 and -2 are chosen arbitrarily. 



Then: 








LG 1 ] = 


[r] [g ][r] T = 


’ 3 


-2 






-2 


4 


[L 1 " 1 ] = 


[r ] [d ] [r] T = 


’ 2 


-2 



-2 4 



[C'] = tr ][C ][r ] T 



4 -4 

-4 4 



The new network is shown in Fig. l-2b. The input impedance is un- 



changed. 



1/2 




1/2 
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It can be seen from example 1-2 that the new network has one more 
element than the original. The circuit configuration and type of elements 
connected between corresponding nodes is also different. 

Guillemin [ 3 ] has shown that an equivalent network can be pro- 
duced by multiplying the same rows and columns of the admittance matrix 
by the same constant. He has also shown that by using a technique 
which seems to be similar to linear programming, an equivalent network 
having all positive elements may be produced from a network with a 
negative element. The following example will illustrate the above 
points . 

EXAMPLE 1-3 

Consider the resistive network of Fig. l-3a. 




The network admittance matrix is 



[G ] = 



60 

-20 

0 



-20 

25 

-10 



0 

-10 

20 



the input impedance Z = ^ ~ ohms . 

If we multiply the second row and second column by 6/5 the new ad- 
mittance matrix is: 



60 


-24 


0 


-24 


36 


-12 


0 


-12 


20 
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The new circuit is shown in Fig. l-3b. 




Fig. l-3b 

Z. = — 7 — ohms 
in 40 

The following example will illustrate the Invariance of transfer 
impedance. 

EXAMPLE 1-4 

Consider the network of Fig. l-4a, connected to a current source. 




The resistances are in ohms, capacitances in farads and inductances 



in henries . 





"0 


0 


o' 




0 


0 


o' 




" 2 


-1 


o' 


[G ] = 


0 


2 


-1 


[C ] = 


0 


1 


0 


[ L _1 ] = 


-1 


1 


0 




0 


-1 


2 _ 




p 


0 


1 




. 0 


0 


0 



The transformation matrix [f ]is arbitrarily chosen as: 



1 0 0 



[ r 3 = 



0 



0 



3 

4 

0 



0 

1 
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Then: 










_ O 




0 0 o' 




’o 0 o" 




2 -2 0 


[G'] = 


« I -f 


[C'] = 


o 

cr> | <sr 
o 


[ L’~ 1 ] = 


-1 i o 

2 4 




.0 -1 *J 




0 0 1 _ 




. 0 0 0. 



The circuit diagram of the new network is shown in Fig. l-4b. The 
transfer impedance has remained invariant. 



2/3 2/3 




I . 3 Realizability 

It has been shown by Guillemin [ 3 ]that if the transformation 
matrix [ F ] is real and non-singular, the positive definite character 
of the stored energy functions and loss function of the network being 
transformed remains unchanged." Consequently, under the condition 
that [F] is real and non-singular, a passive network under the trans- 
formation of (l-3)is always transformed into another passive network. 

J. D. Schoeffler [5 ] has demonstrated that the transformed network 
can always be realized by imbedding the original network in an array 
of interconnected ideal transformers. However, the realization of 
transformed network as a transformerless network with positive elements 
only presents particular problems and imposes some restrictions on the 
values of the elements of the transformation matrix [ F ] . These re- 
strictions and allowable ranges for the values of the elements of [F] 
such that the transformed network can be realized with positive R, L, 

C elements and no mutual inductance is discussed later. 

After the transformation is performed and the admittance matrix of 
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the new network is obtained, it is possible to determine whether the 
new network can be realized by positive R, L, C elements. The follow- 
ing are the necessary and sufficient conditions for [ Y' ] to be realiz- 
able with positive R, L, C elements , Howitt [ 2 ] . 

(1) The off diagonal elements of each of the matrices [G 1 ] , 

[ C' ] and [ L* 3 ] has to be less than or equal to zero, i .e . : 



g 'iJ 




0 


for all i / j 


(l-7a) 


G 'i) 


$ 


0 


for all i / ] 


( 1 — 7b) 


r:‘ 

1J 




0 


for all i ^ j 


(l-7c) 



(2) The sum of all the rows of each of the matrices [G ] , 

[ C ] a nd [ L ^ ] ha s to be non- negative , i . e . : 
n 

L g\. ;> 0 (l-8a) 

3 = 1 13 

n 

I G ' . . >: 0 (l-8b) 

j=l 13 

n -1 

I L\ ^ 0 (l-8c) 

3=1 13 

Conditions (1-7) and (1-8) can be checked by a computer after the 
transformation is performed. 

1 . 4 Transformation of the Branch Matrices 

Since the elements of a network are explicitly presented in the 
branch admittance or impedance matrices, it is sometimes advantageous 
to perform the transformation on the branch matrices . 

The branch current and voltage matrices of a network can be 
written as: 





'•bi 






I b = 


'hi 

• 


II 

XI 

> 


v bi 




l bn 




.V 
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where .^..and ^represent the current through and the voltage across 
branch i respectively. The equilibrium conditions of the network can 
then be written as 

[V b ] = [Z b ]tl b ] (1_6a) 



or 



[i b ] = [v b ][v b ] 



(1 — 6b) 



where [Z^] is the branch impedance matrix formed as follows: 

Each diagonal element Z.. is the total impedance in branch i. 

Each off diagonal element Z, . is the mutual inductance common to 

U 

branches i and j and 



1X3 = Cz. ] 

b b 



-1 



( 1 - 7 ) 



If the loop current matrix of the network is defined by 

h 

tn- 



m 

where I. is the current in loop i then [I 3 and [ I. ] are related by 
i b 



[I b 3 = Card 3 

where [O' 3 is the network tie set matrix formed as follows 

^+1 



(1-8) 



Each element a.. = 
U 



if I,, and I. are in the same direction 
bi i 



if I, , and I. are in the opposite direction 
bi i 



-1 

0 if branch i is not included in loop i 
Similarly if the cut set voltage matrix [ V 3 of a network is de- 
fined by 

' V 



[V 3 = 



1 



V. 



V, 



15 



where V is a voltage with an arbitrary direction assigned to the branch i 
of a tree of a network , then 

[v b := [£ ] T tv ] d-9) 

where [ jS ] is the network cut set matrix formed as follows: 

1 +1 if V, . and V, are in the same direction 
bi 1 

-1 if V^. and V. are in the opposite direction 

0 if branch i is another branch of the tree 
which includes V. 

l 

Schoeffler [ 4 ]has shown that if a matrix [D ] is formed such that 

[n] = [p] T [r][o] (l-io) 

where [ T ] is the transformation matrix defined in section 1.2, then 

instead of applying [T] to the network admittance matrix, we can 

apply [ II ] to the network branch admittance matrix [ Y, ] . 

b 

The branch admittance matrix [Y 1 ^] of the equivalent network is 
then given by 

[Y- b ] = [n][Y b ][n] T (i-ii) 

The advantage of this method is that the elements of the net- 
works are readily seen from [Y' b ] and the need for transformers in 
the synthesis of the equivalent network is apparent from presence of 
off diagonal elements in [Z 1 ^] . The disadvantage lies in the fact that 
[ n ] has to be formed from equation (1-10) which means that cut set and 
tie set matrices have to be constructed before [ H ] can be formed. 

1 . 5 Continuously Equivalent Networks 

J.D. Schoeffler [ 5 ] has developed a method whereby a set of 
equivalent networks is produced as a function of an independent 
variable x. The elements of each member of the set of equivalent net- 
works are functions of the independent variable x. For any value of this 
variable a network equivalent to the original network is obtained. 

Consider that the original [Y ] matrix is Y . Define a new matrix 

o 

Y(x) whose elements are some functions of x such that 

Y(0) = Y 
o 
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( 1 - 12 ) 



Choose r such that: (where [ U ] is the unit matrix) 

in = [u 3 + [B ]Ax 

now 

[ Y' (x) ] = [F][Y(x)][r] T 

= ([U ]+ [B ] Ax) [Y(x)] ([U ] + [B ]Ax) T 

= [ Y(x) ] + ([ Y(x) ] [ B ] T + [ B ] [ Y(x) ] ) Ax 

+ [B ] [ Y (x) ] [ B ] T Ax 2 

If Ax is small then[Y' (x)]will differ.from[Y(x)]by a small amount. 
Define: 

[AY(x)3 = [Y 1 (x) ] - [Y(x)3 (1-13) 

then Equation (1-10) becomes: 

[AY(x)] 

—7 = C B ] [ Y(x) 3 + [ Y(x) 3 [ B 3 + C B 3 C Y(x) 3 [ B 3 Ax 

as Ax 0 

CdY(x)3 

— £ = [B 3 [Y(x)3+ [Y(x)3[B 3 (1-14) 

If the elements of B are chosen arbitrarily^ Equation (1-11) can be 
integrated to yield the elements of[Y(x)3as a function of x. Schoeffler 
integrates Equation (1-11) numerically and changes the assumed values 
of the elements of B as the integration proceeds. In the numerical in- 
tegration process the values of the elements of B are chosen so that 
the new matrix Y'(x) is dominant.* Thus the elements of B are actually 
a point function of x as the integration process proceeds. 

For a network composed of R, L and C we can write 

[Y(x)3 = [G(x)3+ S[C(x)3+ | [L _1 (x)3 (1-15) 

★ 

i.e. the realizability conditions of Sec. 1.3 are satisfied. 
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substituting Equation (1-12) in (1-11) then expanding and equating 
terms of equal powers of S on both sides we obtain: 



d[ ^ x)] = [B 3 [ G (x) 3 + [ G(x) 3 C B ] T 


(l-16a) 


d[ ^ x W:i = [B ][C(x) ] + [ C(x) 3 [ B ] T 


(1-1 6b) 


[L '5x~ 1X)] - [B ) U l « ] + CL ‘m Kb 3 t 


(l-16c) 


which can be shown to reduce to the form 

dC ® (x) 3 = CM 3 [ G (x) 3 


(l-17a) 


= [M 3 [ c (x) 3 


(1 — 1 7 b) 


~d~x^~ = [M ][L " 1(x) ] 


(l-17c) 


or 

- [ M ] [ Y(x) 3 


(1-18) 


where [ M 3 is a matrix whose elements are some combination of the 
elements of [ B 3 . Matrix [ M 3 can be found by solving the matrix 
equations: 

[B 3 C G (x) 3 + L G(x) 3 [ B 3 T = [M 3 L G(x) 3 (1-19) 

or 

[B 3 [ Y(x) 3 + [ Y(x) 3 C B 3 T = [M 3 L Y(x) 3 (1-20) 



i .e . 

[M ] = [B ] + [ Y(x) 3 [ B ] T [ Y~ 1 (x) ] (1-21) 

Substituting Y(0) = Y q in Equation (1-21) yields matrix [M3 which 
can then be substituted in Equation (1-18). After integration this 
yields the admittance matrix [Y(x) 3 of the equivalent networks as 
function of the independent variable x. 
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The Diagonal Variable Transformation and Equivalent 
Networks with Constraints 



2 . 1 Introduction 

In generating equivalent networks using Equation( l-S^or in gen- 
erating continuously equivalent networks by integrating Equation (l-15, ) 
the designer has no criterion for choosing the elements of the transfor- 
mation matrix [ P ] to achieve a design specification such as changing 
the value of an element connected between two nodes in a given manner 
or eliminating a given element from the network. The designer is not 
even able to choose the elements of [rl such that the transformed 
network can be realized by positive R, L, C elements with no mutual in- 
ductance. In this chapter a simplified transformation matrix, namely a 
diagonal matrix, is chosen for the transformation. The acceptable range 
of elements of [r] such that the transformed network can be realized 

by positive R, L, C elements only is found. By choosing the elements 

X x 

of [r] to be exponential functions of a real variable x, i.e. e , a 
group of equivalent networks are generated and it is shown that by 
choosing the values of X 's certain constraints can be imposed on the 
network. However as we shall see in Chapter 4, the diagonal form can't 
always generate equivalent networks. 

2 .2 The Diagonal Transforma tion and Realizability Conditions 

The transformation matrix is considered to be a diagonal matrix with 
positive elements, i.e. 



[r]= [r] T = 




0 

0 



(2-1) 



0 



0 
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where 



y. 

u 



for i = 1 



n 



( 2 - 2 ) 



The admittance matrix of the transformed network then becomes: 

[Y 1 ] = [r][Y ][T] (2-3) 

In general 

i _ i 

(2-4) 



[Y ] = [G ] + S[C ] + | [L X ] 



Substituting (2-4) in (2e3) we obtain 

[Y i ]=[r][G][r]+s[r][c][r] + ^[r][L" 1 ][r] 

(2-5) 

Since in general [ Y' ]will also be composed of R, Land C then 
we can write 



[ Y 1 ] = [ G' ] + S[C’]+ | [ L' 1 '] 

Comparing (2-6) and (2-5) we can write: 

[ G' ] = [T] [G ][T] 

[c 1 ] = [ r] [ c 3 L r 3 
[ l ' -1 3 = [r,][L _ 1 ][r] 



(2-6) 

(2-7a) 
(2 — 7 b) 
(2-7c) 



A typical term of transformed matrices [G 1 ] , [C 1 ] and [ L 1 ^]is 
of the form 



g . = y . . y. . g . . 
u ii jj u 



(2-8) 



By assumption of Equation(2-2)all the y 's are positive. Hence 
the sign of the elements of the transformed matrices [G‘ ] , [ C']and 
[ L' ^ ] remain unchanged. If the original network was made up of 
passive elements the off diagonal elements of the [G ] , [C ] and 
[L ■*■] matrices were negative. Hence the off diagonal elements of 
the transformed matrices are also negative and one of the realizability 
conditions stated in section 1.4 namely Equation (l-6)is always sat- 
isfied and the circuit elements connected between any two nodes other 

than node to datum, remain positive after the transformation. 

th 

The sum of the i row of the transformed matrices [G 1 ] , [ C ' ] 
and [ L' are given by the following equations: 
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n 


n 






S'.(G) = Eg 1 ,, 

1 j=l 13 


= E 
j=l 


y. . y. . g . . 


(2-9a) 


n 


n 






S'.(C) = SC'.. 

1 j=l 13 


= E 


y.. y.. c.. 
n Ji n 


(2 — 9b) 


j=l 




s , .(l" 1 )= el" 1 

1 j=l 13 


n 


T -l 
y. . y. . L, . 
n '}) i) 




= E 
1=1 


(2-9c) 



It is clear from the method of the construction of the admittance 

matrix that the above sum of the rows corresponds to the value of the 

element connected between node i and datum, i.e. g'. , C'. and 

10 10 

L' . 

10 

If the second of the realizability conditions that is Equation (1-7) 
is to be satisfied we must have: 

n 

S' . (G) = I y.. y.. g.. ^ 0 for all i (2-10a) 

1 j=l 11 33 13 

n 

S'.(C) = I y. . y. . C, . > 0 for all i (2-10b) 

1 j=l 11 J 1 1J 

-1 n -1 

S'.(L ) = E y. . y. .L.. ^ 0 for all i (2 -10c) 

1 11 'n iJ 

Equations (2 -1C) represent a set of constraints on the elements of the 
transformation matrix. Since^from Equation (2-2) y^ is always positive 
it can be cancelled from the set of Equations (2-10). Thus in general 
there are 3n linear equations that have to be satisfied, n equations for 
each of the [ G 1 ], [C 1 ] and [L 1 matrices comprising [Y 1 ]. If 
Equations (2-10) are satisfied with the equal sign then the corresponding 
element connected between the corresponding node and datum is zero. 
For example, if S'^G) = g^ is zero then the conductance connected 
between node 1 and datum is eliminated from the circuit. From the 
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above we can conclude that there are no set of y 's which would satisfy 
all of the 3n equations of (2-10) with the equality sign, since this 
will result in an equivalent network with no node to datum element^ , 
i.e. with infinite input admittance at any port. 

The 3n inequalities of the form (2-10) together with the n in- 
equalities of the form (2-2) are a set of 4n inequalities which must be 
satisfied by any set of y's to assure the fact that the transformed net- 
work is realizable with passive elements only. Obviously these in- 
equalities are all satisfied if all of the y 's are unity, since the original 
networks results in this case. 

From the discussion on Chapter 1 it is obvious that to keep the in- 
put imittance between node k and datum constant y must be equal to 

K. 

unity and to keep the transfer impedance between parts k and l con- 
stant we must choose y and y to be unity. There are then n-1 or n-2 

K JC 

remaining y 's to be chosen satisfying Equations (2-2) and (2-10). For 
example in the case of transfer function equivalence there are a total 
of (n-2) + 3n = 4n - 2 inequalities to be satisfied by n-2 unknown 
y 's. These inequalities are considered in the next section to determine 
if a solution is possible, and if so to ascertain specific numerical ranges 
for the y 's over which the inequalities are satisfied. These ranges 
define the domain of the set of y 's for which the transformation of (2-1) 
generates networks with passive components for which the transfer 
impedance between ports i and j remains invariant. Obviously a 
solution cannot always be expected since the transformed network is 
confined to the class obtainable by the diagonal transformation matrix 
(2-1). The general case where a complete n x n transformation matrix 
is used will be discussed in Chapter 3. 

2 . 3 The Domain of the Solutions 

Each of 3n inequalities of Equation (2-10) represents a linear 
constraint on (n-1) or (n-2) y 's depending on whether the input im- 
pedance or transfer impedance is to remain invariant under the trans- 
formation. When any of the inequalities (2-10) are satisfied with the 
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equality sign the resulting equation represents a hyperplane in (n-2) 

dimensional y - space.* This hyperplane divides the (n-2) dimensional 

y - space into two distinct regions. The points on one side of the 

hyperplane satisfy the corresponding inequality. For example, if the 

inequality g'. > 0 is used the equation of the hyperplane is g' =0. 

lq lO 

Points on one side of the plane yield g'^ Q greater than zero and points on 

the other side yield g 1 . less than zero.When all of the inequalities 

10 

(2-10) are considered-, ‘the domain of the allowable set of y 's is de- 
termined by a hypervolume in y-space bounded by a set of at most 
(4n-2) hyperplanes , 3n of which are of the form; 

n 



Z y. . g . . = 0 
_ n JJ iJ 



( 2 - 11 ) 



and the remaining number in the form 

y.. = 0 (2-12) 

ui 

where i ^ k and i ^ l since y, . = y =1 due to the required in- 

KK <vJo 

variance of the transfer function. Points on the planes representing 
Equation (2 -12) are not acceptable since these points violate the required 
non-singularity of the transformation matrix (2-1). 

Any point on the hypervolume described above represents a set of 
y 's which would satisfy all of the inequalities of (2-2) and (2-10). 

The realizability volume may be determined in a systematic manner 
as follows. Solve all the possible combinations of the (n-2) simultaneous 
equations, taken from the set of (4n-2) equations of the form (2-11) and 
(2-12) for a set of (n-2) y 's, and check whether these values also 
satisfy the remaining (4n-2) - (n-2) = 3n inequalities which are not 
used in choosing the (n-2) simultaneous equations. Any set of y 's 
which satisfy these conditions is a vertex of the realizability volume. 

The hypervolume is bounded on the sides by any of the hyper surfaces 

* 

For convenience of writing, the case of transfer impedance equi- 
valence is considered. The case of input impedance equivalence follows 
in the same manner. 
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which passes through (n-2) vertices. 



In general a total of 



(4n - 2) ! 



(4n - 2)! 



(2-13) 



[ (4n - 2) - (n - 2)] ! (n - 2) ! 3n!(n-2)! 

linear simultaneous equations must be solved and checked in the above 
manner. In practice many of the solutions are trivial and the process 
may be programmed on a digital computer. After all vertices of the 
realizability hyper-volume are determined, the boundary hyper-planes may 
be determined from the simultaneous equations whose intersections 
determined the vertices. The point where all of the v 's are equal to 
one always satisfies the realizability conditions and is always in- 
cluded in the realizability hyper-volume. 

If there is no set of y 's which would satisfy the Equations (2-2) 
and (2-10) we can deduce that the transformation of (2-1) will not pro- 
duce any equivalent network with passive components, and the real- 
izability hyper-volume in this case is reduced to a single point where 
all of the y 's are unity. If a realizability volume does exist it is possible 
to remove at least one node to datum element by choosing a set of y 's 
represented by a point on one of the boundary hyper-planes of the form 
(2-11). 

EXAMPLE 2.1 



Consider the RLC network shown in Fig. 2.1, 




where all the conductances are in mhos, all the capacitances in farads 
and all the inductances in henries . 
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-1 


0 


0 ' 






-1 


1 


0 


0 




[G ] = 


0 


0 


1 


0 


[C ] = 




0 


0 


0 


1 





0 0 0 

1 -1 0 

-13 0 

0 0 2 



C L _1 ] = 



0 0 

0 1 

0 -1 

0 0 



0 0 

-1 0 

4 -1 

-1 1 



It is desired to keep the transfer impedance invariant so let: 





’ 1 


0 


0 


o 


1 — 1 

TJ 

l 1 

II 


0 


y 22 


0 


0 




0 


0 


CO 


0 




o 


0 


0 


1 



The transformed matrices [G'] , 



[ C' ] and [ L' ^ ] are: 



2 



[G'] = 





0 



0 0 
0 0 




0 



0 0 



0 1 



2 



0 0 0 



0 



[C 1 




0 



y 22 V 22 y 33 

~ y 22 y 33 3y 33 

0 0 



0 

0 

2 
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[L -1 ] = 



0 

0 

0 

0 



y 22 

" y 22 V 33 



" y 22 y 33 



4 y 



33 



-y. 



33 



-y. 



33 



The inequalities (2-10) in this case are 12 inequalities but only 8 of 
them contain the variables y and y and of these only 6 are distinct and 

6 O 

are tabulated below . 

From [ G ’ ] we have 



(a) 2 - y, 
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:> 0 

(b) -1 + y 22 £ 0 

From [C']we have 

(c) y 22 -y 33 » 0 

(d) -y 22 + 3y 33 i 0 

From [ L' ^ ] we have 

(e) -y 22 + 4y 33 - 1 



^ ~ y 33 + 1 ^ ° 

In the y-plane the above inequalities define regions as shown in Fig. 2.2 
Taking into account the Equation(2-2)which states 



>22 > ° 



>33 > 0 



the realizability hypervolume in this two dimensional example will be a 
four sided polygon shown shaded in Fig. 2.2, 

To solve the above problem with a digital computer we must solve 
all the combinations of the six equations obtained from the inequalities 
a , b , c , d , e and f two at a time . 

Table 2.1 presents all the possible solutions of the group of two 
simultaneous equations together with the letter representing the solutions 
on the y-plane of Fig. 2.2. The solutions which yield valid inter- 
sections which satisfy all of the other inequalities are indicated by the 
asterisks and are the vertices of the realizability area. The boundaries 
of the volume are given by the lines representing a, b, e and f. 
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TABLE 2 . 1 



EQUATIONS 




SOLUTIONS 


INTERSECTION 


PAIR 


y 22 


y 33 


POINT 


a b 


none 


none 


none 


a c 


2 


2 


A 


a d 


2 


2/3 


B 


a e* 


2 


3/4 


c* 


a f* 


2 


1 


D* 


b c* 


1 


1 


E* 


b d 


1 


1/3 


F 


b e* 


1 


1/2 


G* 


b f* 


1 


1 


E* 


c d 


0 


0 


0 


c e 


1/3 


1/3 


H 


c f* 


1 


1 


E* 


d e 


3 


1 


K 


d f 


3 


1 


K 


e f 


3 


1 


K 
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Fig. 2.2 'Realizability regions in the % plane. 



2 .4 The Exponential Searching Function 

If a volume of realizability does exist then for each point in this 
volume there exists an equivalent network. These equivalent networks 
can conveniently be investigated by letting the elements of our trans- 
formation matrix [ T ] to be exponential functions of an independent 
real variable x. Thus let 

y = e \i x (2-14) 

li 

where X„ is a constant which may be chosen arbitrary or in accordance 
with some additional constraint placed upon the desired network. 

The exponential function of Equation (2-14) satisfies the require- 
ment of Equation (2-2) in that it is greater than zero for all finite 
values of the independent real variable x. Hence the transformation 
matrix [P] remains real and non-singular for all finite values of x. 
Moreover when x = 0, y = 1 , and for this value of the variable x the 
original network is reproduced. 

When the exponential functions of Equation (2-14) are used as 
elements of the transformation matrix, [T] becomes 



e X ll X 0 . . . 0 

0 e X 22 x 0 



[ r (x) ] = 



(2-15) 



.^nn x 



A typical element connected between two nodes in the transformed 
network becomes: 



g 



ij 



% 



X. .)x 
JJ 



U 



(2-16) 
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And a typical node to datum element becomes: 



n 

g'. = I g.. 

10 j=0 1J 



= L g. . e 

J-l 1J 



( X.. + X. ,)x 
11 JJ 



(2-17) 



And a typical equation of the hyperplane used in determining the real- 
izability volume becomes: 



n X. . x 

L g. . e 11 =0 (2-18) 



g'.. and g'. can easily be plotted by a digital computer as a function 
ij 10 

of x. The range of x for which all of the 3n equations of the form of 
Equation (2-17) remain non zero is the range for which transformation 
of (2-15) yields equivalent networks with positive R, L, C elements. 
Differentiating (2-15) with respect to x yields: 



where: 



d [r(x)3 

dx 



[r] [x ] = [x ] [r] 




0 




0 

0 



[x] = 



, 0 



0 



# # 



• # 



X 



nn 



Differentiating (2-3) with respect to x yields: 



d[ Y'] 
dx 



d[T] 

dx 



[y ][ r]+ [r] [ y ] 



d[r] 

dx 



(2-19) 



( 2 - 20 ) 



( 2 - 21 ) 
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Substituting (2-19) in (2-21) we obtain: 



d[Y‘] 

dx 



[ X ] [ Y' ] + [ Y‘ ] [ X ] 



( 2 - 22 ) 



Thus the use of the exponential search function for generating equivalent 
networks is the same as integrating Equation (2-22) with the initial 
condition that [ Y' ] = [ Y ] when x = 0 . 

Comparing Equation (2-22) with Equation (1-11) we observe that 
the two are identical with [ B ] replaced by [ X ] . 

The choice of the X 's may be made to impose additional constraints 
upon the transformed network as follows . 

(a) To keep the transfer impedance between port k and the output 
voltage between nodes A and m constant (no common node 
between input and output) set X, , = X 00 = X =0. 

(b) To keep the admittance between nodes i and j invariant under 



the transformation set X.. = -X.. 

n JJ 

can be seen that in this case 



(c) From Equation (2-16) 



. From Equation (2-16) it 

g\ . = g. . . 

ij ij 



g'.. = 

ij 



g.. e 

IJ 



( X.. + X..) X 
ii JJ 



(2-23) 



and 



(X 



kje 



= g k* e 



kk 



+ V X 



i 

(2-24) 



Dividing (2-23) to (2-24) we obtain 



U 



/ ■ 



9 u 



/ g. 



[ ( x.. + x..) - ( x, , + x 

n jj kk 






Hence to keep the ratio of the admittance connected between 

nodes i and j to the admittance connected between nodes k 

and A constant set ( X. . + X. .) = ( X, , + X.J . 

ii jj kk AA 
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(d) Since the circuit elements between nodes i and j increase or 

decrease monotonically with x depending upon the sum 

(X.. + X.. ), the X 's can be chosen so that a specified com- 
11 JJ 

ponent changes more rapidly or slowly compared with others 
over the range of x. 

The use of the exponential search function and constraint (b) for choosing 
the X's is demonstrated in example 2.2. 

EXAMPLE 2.2 




It is desired to keep the transfer impedance as well as the admittance 
connected between nodes 2 and 3 constant. 

From Fig. 2.2 it can be seen that since x = 0 corresponds to 
point E, to stay in the area of realizability we have to choose and 
Xgg such that 7^2 increases andy^ decreases with x. Hence X ^ 
must be negative and X 22 m ust be positive, i.e. X 22 = "^33 = constant 
the value of the constant depends on the desired rate of change of the 
circuit elements . Since at the moment we are not concerned with rate 
of change of elements we choose 



y 22 y 33 1 



The [G'l , [C']and [ L' matrices would then become: 
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L G‘ ] = 



[C'] = 



-e 
0 
0 

2 
0 
0 
0 

0 
0 
0 
0 

The new circuit elements become 

g 



[L'] = 



X 

-e 


0 


2x 

e 


0 


0 


-2x 

e 


0 


0 


0 


0 


2x 

e 


-1 


-1 


3e -2x 


0 


0 


0 


0 


2x 

e 


-1 


-1 


/I " 2x 

4e 


0 


-x 

-e 



0 

0 

0 

1 

0 

0 

0 

2 

0 

0 

-e 



10 



12 



20 



y 30 
9 '40 “ 1 



2 -e X 


c 


X 

e 


c 


2x x 

e - e 


c 


-2x 

e 


c 



c' 



10 

20 

23 

30 

40 



= 2 



2x , 
e - 1 



= 1 



3e- 2x - 1 



= 2 



T ' = 

L 20 



l/(e 2X - 1) 



L '23 = 1 

L’ 30 = l/(4e" 2x -e‘ x - 1) 
L 34 = 6 

l ' 40 = 1/d - e_X ) 
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The variation of the conductances, capacitances and the inverse in- 
ductances with the independent variable x are shown in Figs. 2.4, 2.5 
and 2.6 respectively. 
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X 



Fig. 2.4 Variation of Conductances with x, Example 2.1 
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The General Transformation 



3 . 1 Introduction 

In Chapter 2 a diagonal transformation matrix was used and it was 
shown that the realizability hyper- volume , when it exists, can be 
determined by the manner described in that chapter. The diagonal trans- 
formation matrix however is limited in application. For example it is 
obvious from Equation (2-8) that if no conductance existed between nodes 
i and j in the original network, none of the equivalent networks gen- 
erated will have a conductance between nodes i and j. But the sim- 
plicity of the transformation (2-3), linearity of the realizability 
equations (2-10), and the ease of applying constraints on the transform- 
ed networks warranted the use of the diagonal transformation. 

In this chapter the full n x n transformation matrix [T] is used. 

The equations for the elements of the transformed networks are obtained. 
The realizability conditions and the allowable range for the elements of 
the transformation matrix [T] so as to yield to equivalent networks 
with positive R, L, C circuit elements is then discussed. 

3 . 2 The General Transformation and Realizability Conditions 

The transformation matrix r is chosen to be 










(3-1) 




y 



nn 



A typical element of the conductance, capacitance, and inverse 
inductance matrices [G 1 ], [C‘] and [ L 1 of the equivalent networks 
produced will then be of the form: 
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(3-2) 



n 



g’ = £ 



u 



je,k=i 



y u y jk g i k 



Since the original network is made up of bilateral elements, the 
[G ] , [C ] and [L ^ ] matrices are symmetric. Hence Equation (3-2) 
can be written as: 



n 



n 



g'.. = £ g. . y .. y .. + £ g [y. y. +y, y. ] 

lj k=1 kk lk jk p=1 y pq ' ip jq iq ' jp 

q=p+l 



(3-3) 



Where i = j, Equations (3-3) yield the diagonal elements of the matrices 
[G‘ ] , [ C' ] and [L'" 1 ] . Thus: 



n 2 

g'.. = £ g y + 2 

n , , kk ik 

k=l 



n 

£ 

P=1 

q=p+l 



g y . y . 
y pq ip iq 



(3-4) 



We recall from Chapter 1 that realizability with passive components 
required that the off diagonal elements of the [G 1 ] , [C 1 ] and [L 1 ^ ] 
be zero or negative and that the sum of terms in any row of the above 
matrices be zero or positive. Thus realizability with passive elements 
in the case of the full transformation matrix [T ] requires that 



g'.. ^ 0 

u 

n 

g'. = £ g’.. 

10 j=l 1J 



for i ^ j 



(3r5) 

(3-6) 



For an n port network each of the matrices [ G‘ ] , [ C' ] and 

[ L' is a symmetrical n x n matrix. Hence we have in general 

3n(n-l)/2 inequalities of the form (3-5) and 3n inequalities of the form 

(3-6) to be satisfied. Since the transformation matrix [T] need not be 

2 

symmetrical it contains n elements. The invariance of the input or 
transfer impedance requires that one or two columns of the transformation 
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matrix [ T ] be unit vectors (i.e. zeros everywhere except on the diag- 
onal position). Hence if the invariance of the transfer impedance is 

required a total of 3n(n+l)/2 inequalities have to be satisfied by 
2 

(n -2n) = n(n-2) y 's if the transformed network is to be realizable with 
passive components. 

Examination of (3-3), (3-4), (3-5) and (3-6) shows that the in- 
equalities (3-5) and (3-6) which have to be satisfied are non-linear 
functions of the y 's. Hence the method used in Chapter 2 for the diag- 
onal transformation matrix is not applicable here, since the equations to 
be solved are not simultaneous linear equations. To be able to apply 
the method of Chapter 2 to the transformation of (3-1) used in this chap- 
ter, we will linearize the inequalities of (3-5) and (3-6) as shown in 
the next section. It must be noted that in the case of transformation 
with the full transformation matrix of (3-l),in contrast to the case of 
diagonal transformation matrix, the y 's can take both positive and neg- 
ative values . 

3.3 Linearization 

We know from the discussion in Chapter 2 that the inequalities 
(3-5) and (3-6) are satisfied when y_ = 1 and y = 0 for i / j . 

' This will result in the transformation matrix [ F ] to be the 
identity matrix [ I ] and the transformation will yield the original net- 
work . 

We can linearize the inequalities (3-5) and (3-6) by expanding the 

matrix [Y 1 ] in a Taylor series expansion about the point y = 1 and 

y..= 0 i 4 j as shown below: 
ij 

From (1-3) we have 



[ y 1 ] = [r] [y ] [r] T 


(3-7) 


or from (l-6a) we have 




[G 1 ] = [r] [g ] [r] T 


(3-8) 
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Hence 



[AG 1 ] = [Ar][G ][r] T + [r][G ][Ar] T 



n 



[AG’] = [ E 



i,j=l 7 ij 



I to HI]+ [r][G ] 

sy,. u 



n 

x / I 



Ay.. 

l.J-1 av ii U 



(3-9) 



But 



sen 

6r ij 



4? ii " 



Ay. . 
ij 



Hence 



n 

E 

i =1 
j=l 



or 

dy 



ij 



Ay . . 
U 



Ay ll 


Ay 12 


• 


• 


• 


Ay m 


Ay 21 


Ay 22 


• 


• 


• 


• 



Ay 



nl 



Ay 



nn 



=[Ay ] 



(3-10) 



Since we are expanding [G'] with a Taylor series about the point 

y . , = 1 ,y..=0 i^i we can write 
n ij 



v . . = 1 + 6 . . 

ii ii 

y . . = 0 + 6 . , 
'ij U 



(3-11) 

(3-12) 
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i.e. 



where 6.. and 6. . are incremental change on y , . and y, , , 
n ij ii ij 

6 .. = Ay., and 6 .. = Ay . and the transformation of (3-1) with 
n ii ij ij 

6.. = 6.. = 0 yields the original network. The origin of the 6*space 
corresponds to [ T ] being the identity matrix. 

Substituting (3-11) and (3-12) in 3-9 we can write 

[6] = [Ay] = I ~ Ay.. (3-13) 

Substituting (3-13) in (3-9) we obtain 

[AG'] = [ 6 ] [ G ][r] T + [ T ] [ G ][6] T (3-14) 

Expanding [G‘] with the Taylor series around [G ] yields 

[G 1 ] - [G'] | q + [AG 1 ] | q (3-15) 

where [G 1 ] I and [AG 1 ] I denotes the matrices of [G']and 
'o o 

[ AG' ] where y. . = 1 and y . . = 0 i^j or where 6 .. = 6 . . = 0 . 

n ij n ij 

From (3-8) 

[G 1 ] | Q = [G ] (3-16) 

Since [r]| = [r] T | = [U] where [U] denotes the identity 

o o 

matrix then from (3-14) we have 



[AG’] | o = [6 ] [G ] + [G ] [6 ] T (3-17) 

Substituting (3-17) and (3-16) in (3-15) yields 



[G'] = [ G ] + [6 ][G ] + [ G ] [ 6 ] T (3-18) 

Hence a typical term of the [G 1 ] , [C 1 ] and [ L' 1 ] matrices would 
be 



n 





+ I 
k=l 



( 9 ,.. 6 



kj ik 



g. 



ik 






where only the first two terms of the Taylor series are taken. 



(3-19) 
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Equation (3-19) will enable us to linearize the inequalities of 
(3-5) and (3-6) and hence apply the method of Chapter 2 to obtain the 
area of realizability and the range of y 's for which the transformation 
of (3-1) yields equivalent networks with passive elements. In essence 
this procedure identifies a small region in 6- space, any point of which 
corresponds to a set of y 's for which the transformation of(3-l)may be 
applied to yield equivalent networks with passive elements. Once 
this region is defined a searching function such as the exponential 
searching function of Chapter 2 can be used to generate equivalent 
networks whose circuit elements are functions of an independent variable 
such as x and the complete inequalities of (3-5) and (3-6) investigated 
to yield equivalent networks with passive elements. 

3 .4 The Exponential Searching Function 

We use an exponential searching function for the elements of the 
transformation matrix [ I ] in the following form: 

X.. x 

y = e 11 (3-20) 

X.. x 

y.. = 1 - e 1J i ^ j (3-21) 

Approximating (3-20) and (3-21) by the first two terms of their Taylor 
series expansion and comparing with Equations (3-11) and (3-12) respec- 
tively we obtain 



X.. x = 6.. (3-22) 

n li 



-X.. x = 6.. (3-23) 

i] iJ 

Equations (3-22) and (3-23) in conjunction with the linearized form 
of Equations (3-5) and (3-6) can be used to choose values of the X 's 
and the range of x, as best illustrated by the following example. 
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EXAMPLE 3.1 



Consider the RC network of Fig .3.1. 




NAANWVv 



(2) 



Fig. 3.1 



It is desired to find the equivalent networks containing passive elements 
which have the same input impedance as the above network. 

The [G ] and [C ] matrices of the above network are: 



[C ] = 



o 

1 — ' 


1 

o 


[G 3 = g 


' 1 


-1 ' 


0 


c 2j 




-1 


1 



Starting with the diagonal transformation matrix 



[r ] = 



i 

o 



22 



The [G'] and [C 1 ] matrices will be 





1 


-y 22 




c i 


0 


[ G 1 ] = g 




2 


lc-3 = 


_0 


_ 2 




_ _y 2 2 


V 22 




C 2 y 2 2 



The realizability conditions are 



1 - y 



22 



?> 0 



-! + y 22 £ o 



C 2 y 22 



^ 0 
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It is obvious that the only value of y 22 which satisfies the above 
equations is y ^ = 1 which means that the transformation matrix is the 
identity matrix. Hence the diagonal form of the transformation matrix 
will not generate any equivalent network. 

Now we use the full transformation matrix 



[r] = 



1 

0 



V 



12 



V 22 



The [ G ' ] and [ C ' ] matrices will be 



C G' ] = g 



V 22 V 12 ' V 22 



V 22 y 12 ' y 22 




C 1 + C 2 y i 2 



C 2 y 12 y 22 



[C'] = 



C 2 y 12 y 22 



^ 2 
C 2 y 22 



The realizability conditions (3-5) and (3-6) become: 

(a) g' 10 = (l-2v 12+yi 2 2+ v 12 y 22 -y 22 )g > 0 



ft >) 9' 20 * y 22 y 12 " y 22 + y 22 * 9 ® 0 
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0 



(c) g, 12 = ( y 22 y !2 - y 22 )g * 



(d) C- 10 = C 1+ C 2 (y 12 y 22+ y 12 ) * 0 



^ C 20 C 2^ y 12 y 22 + y 22 ^ ^ 0 



(f) C' 12 C 2 y i2 y 22 ^ 0 



The linearized form (3-19) of the above equations are 
(a) - 6 12 - 6 2 2 ^ 0 



<b)«i2 + 6 22 * 0 



(c) 1 + 6 22 ‘ 6 12 » 0 



(d) Cj + C 2 6 12 £ 0 



(a) C 2 (6 12 + 2 6 22 + l) J 0 



(f) C 2 6 12 « 0 



Each of the above inequalities when satisfied with equal signs will be 
presented by a line on 6^ 2 , 6 22 space as shown in Fig. 3.2. The 
linearized realizability region reduces to the portion of the line 
6 99 +6 19 =0 in the second quadrant of the 6 space shown on Fig. 3.2 
Using the exponential search function, Equations (3-20)and (3-21) 
become 



v 22 = e 



X 22 x 



y 12 = 1 - e 



X 12 X 
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Equations (3-22) and(3-23)become 



X 22 X 6 22 



' X 12 X 6 12 



From the conditions of realizability for this case we have 5 ^^ = -6 2 = 
constant. The value of the constant determines the rate of change of 
the numerical values of the circuit element. Since we are not concerned 
at this moment about the rate of change of elements we arbitrarily 
choose \^2 = = 1 f° r x >0. 

Hence 



*22 =e 



*12 = 1 ’ e 



X 



Substituting the above values in the expressions for g's and c's yields 



9 10 0 



C 'lO = C l + C 2 (1 - eX) 



9 20 ° 



g'i2 = - ge 



2x 



C ' 2 0 = C 2 e 
c, 12 = c 2 e X (l-e X ) 



From the above equations it can be seen that the range of x for which the 
realizability conditions are satisfied is given by the two equations 

C x + C 2 (l - e X ) £ 0 



Hence 



C 2 e x (l - e x ) <0 



« e x s 1 + -r 1 - 



or 



x in 



C l +C 2 
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The equivalent network is shown in Fig. 3.3 below and the variation of 
the network elements with x are shown in Fig. 3.4 for the case of g = 1 , 



C 2 e (e -1) 
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I 

I 

i 
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Fig. 3.2 Linearized realizability conditions in the 8 plane. 




Fig. 3.4 Variation of Circuit Elements of Example 3.1 
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Conclusion and Suggestions for Further Research 

4 . 1 Conclusion 

In this paper the congruent transformation is applied to the admit- 
tance matrix of an n+1 node common datum network to generate equivalent 
networks which have the same input or transfer admittance. It is shown 
that a transformation matrix whose elements are exponential functions 

XrX 

of a real variable x, such as e , can be used to generate a grouh of 

equivalent networks whose circuit elements are functions of the variable 

x. The range of variation of the elements of the transformation matrix 

required to generate equivalent networks with passive circuit elements 

is found by solving a number of sets of equations called the realizability 

equations. If a diagonal transformation matrix is used these equations 

are linear but when a full transformation matrix is used these equations 

are non-linear, and have to be linearized. If a diagonal transformation 

matrix ^whose elements are exponential functions of the real variable x 

is used, then certain constraints can be placed upon the generated 

equivalent networks by choosing certain relations between the constant 

X x 

factors ( X ' s) in the exponent of the exponential function e 

4 .2 Suggestions for Further Research 

It is believed that the subject of equivalent networks has a great 
potential for application in the field of network synthesis and analysis. 
The following are some ideas for further research on this subject. 

(a) It may be possible to find the general form of a transformation 
which would take a given form of a synthesised network to another form, 
For example it may be possible to find a transformation matrix [T] which 
transforms the Foster form of a synthesised network to Cauer's form for 
the same input impedance . 

(b) The effect of the presence of mutual inductance in a network 
on the transformation can be investigated with the aim of transforming 

a network with mutual inductance to a network without mutual inductance. 

(c) Finally, from curves of Fig. 2.4, 2.5 and 2.6^it can be seen 
that if the values of the components of the network vary according to 
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the curves shown, then the transfer impedance and input impedance 
remain invariant, i.e. the multiparameter sensitivity of the input or 
transfer impedance is then zero. Hence, multiparameter sensitivity 
can be studied from the point of view of equivalent networks. 
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Appendix I 

Invariance of Input Impedance 



It was stated in Chapter 1 that if the column 1 of the transformation 
matrix [ 1 ] is a unit vector (zero everywhere except for a one at the i** 1 
row) then the input impedance at port i, that is Z . , remains invariant 
under the transformation. In this appendix a proof of this fact based on 
matrix portioning is presented. 

Consider an n+1 node network with a current source connected be- 
tween each of the n nodes and the datum node as shown in Fig. 1-1 . 



1 



3 




The equilibrium equation of the network can be written in matrix 
form as 



[I 3= [Y 3 L V ] 



(1-1) 



with [ I ] , [ Y 3 and [ V ] as defined in Chapter 1 . 

Expanding the matrix Equation (1-1) and using Crammer rule, the 
voltage between node i and datum is given by 



A. 



V. = 
i 



il 



A. 



r i + 



i2 



V- • 



A.. 

+ 

A i 



in 



n 



( 1 - 2 ) 



where: 

A = the determined of the admittance matrix [ Y ] 

^ij = m * nor f° rme d from Aby deleting row i and column j 

If the network is excited only at node i by a current source I then* 

i 

I . = 0 for i = 1 . . . . n i ^ j 
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Hence 



(1-3) 




and 

A.. 

Z. = (1-4) 

1 A 

The invariance of the input impedance at port i under the congruent 
transformation 

[Y']= [T] [Y ] [r] T (1-5) 

with column i of [ F ] being a unit vector would be proved if we can 
show that the ratio A„ / A of Equation (1-4) remains invariant under 
transformation (1-5) . 

For simplicity of proof we consider the invariance of the input im- 
pedance at node 1. Since by remembering the nodes we can call any node 
node 1 , there is no loss of generality and the proof would be acceptable 
for the invariance of the input impedance at any port. 

Consider a real non-singular transformation matrix [T] where 
column 1 is a unit vector 



[r] = 



0 y 



12 



22 



• # 



• • 



• y 



In 



. y 



2n 



0 y 



n2 



• • 



. y 



nn 



(1-6) 



The admittance matrix [ Y'] of the transformed network is given by 

[ Y' ] = [1 ][Y ][r] T (1-7) 

where [ Y ] and [ Y‘ ] are as defined in Chapter 1 . 
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Since the determinant of the product of matrices is equal to the 
product of the determinants [P] then: 

A' = | T | 2 A (i-8) 

where A' = determinant of [Y 1 ] 
and | T | = determinant of [ I ] 

We can partlttn the [ T ] , [ T ] T , [Y ], [ Y' ] matrices as shown 



below. * 



wciu w • 

y, ll '• • • • y 'ln 




1 I Y 12 * * ‘ V ln 






-i 


• | 

| X' 

1 

• | 

y' ' 

y nl ! 




0 1 

. 1 a 

( 

I 

• 

f 

0 t 

i 

i i 

u „ 


y ll ! ‘ y ln 




O 

• 

• 

• 

O 

i— i 


__ _ 1 - _ _ - — - - — — ^ 




I 

• 


• 1 

1 

1 X 

. y n ! 




a i2 1 
1 

1 t 

. 1 0 1 

a *„ ! 

\ 

„ \ 



(1-9) 

Performing the matrix operation we obtain 

[ X' ]= [ a ] [ X ][a] T (I- 10) 

Hence the admittance matrix of the new network is given by: 

★ 

y. . i ^ j = the admittance connected between nodes i and j 

ij 

y,. = the sum of all admittances connected to node i 

li 
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[Y'] = 




\ 

I 



'■ Ca][X ][cy] T 



I 

\ 



I 




I 

i 



(1-11) 



The input impedance at node i of the new network can be written 
from Equation (1-4) to be 



7 . = A 'll = det. ( [g] [X ] [g] T ) 

1 A' A' (1-12) 

Since determinant of product of matrices is equal to the product of the 
determinants [ 7 ] we can write • 



det ( c Of 3 C x ][o 3 T ) = I a | 2 I X |= I a | 2 A n (1-13) 

where | a | = determinant of [ a ] . 

Since the 1st column of [ I ] is a unit vector 

| r | = |<y | (1-14) 

Substituting (1-14) in (1-13) we obtain 

det ( [ cy ] C cy ] [a ] T ) = | T | 2 A x 1 (1-15) 

Substituting (1-8) and (1-15) in (1-12) yields 

z , = 11 \ A n ^i 

1_ in 2 A " A 

Hence ‘ 

Z 1 = Z x (1-16) 
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Appendix II 

Computer Program for Example 2.1 



• • JOB0184F » ARDALAN 

PROGRAM EQTRANF 1 

DIMENSIONXi 500) ,010(500 ) ,612(500 > ,G20( 500 ) ,630 < 500 ) ,G40( 500 ) .C10( 
1500) *C20(500) .C23( 500) »C 30(500) ,C40<500) .P20(500) .P23(500) .P30(50 
2 ) .P 34 (500 ) ,P40( 500) . ITITLE(12) »GA2 ( 500 ) .GA3( 500) 

GM2 = 1 • 

GM3*-1 • 



X( 1 )“0. 
A=0 • 8 
DO 10 
G12 ( I 
G10 ( I 
G20 ( I 
G30 ( I 
GAO ( I 
CIO ( I 
C23 ( I 
C20 ( I 
C30 ( I 
C40 ( I 
P23 ( I 
P20 ( I 
P34 ( I 
P30 ( I 
P40( I 
6A2 ( I 
GA3 ( I 
X( 1 + 1 



10 



22 



23 



24 



1=1,500 

“EXPF (GM2*X( I ) ) 

*2 »-G12 ( I ) 

“EXPF ( 2 • *GM2*X ( I ) )-G12( I ) 

“EXPF (2 ,*GM3*X ( I ) ) 

= 1 . 

* 2 • 

=EXPF ( ( GM2+GM3 ) *X ( I ) ) 

=EXPF(2.*GM2*X< I ) )-C23( I ) 

=3,*EXPF( 2.*GM3*X( I ) )-C23( I ) 

*2. 

=EXPF ( ( GM2+GM3 ) *X ( 1 ) ) 

“EXPF ( 2 ,*GM2*X ( I ) )-P23 ( I ) 

“EXPF (GM3*X< I ) ) 

“4**EXPF( 2,*GM3*X( I ) ) -P23 ( I ) ~P34( I) 

“1.-P34 ( I ) 

“EXPF (GM2*X( I ) ) 

“EXPF ( GM3*X ( I ) ) 

=X( I >+( A/500. ) 

CONTINUE 

DO 22 1*1.12 

ITITLEt I ) “12 

I T I TLE( 1 ) “8HARDALAN 

LABLE“4HG10 

CALL DRAW (500.X.G10.1 .0. LABLE , ITITLE.0.0,0.0.0.0.5.5 ,0. LAST ) 

DO 23 1=1,12 

I T I TLE( I ) =12 

I T I TLE ( 1 ) “8HARDALAN 

LABLE=4HG12 

CALL DRAW( 500.X.G12.2 .0, LABLE, I T I TLE .0,0 .0,0 . 0 .0. 5. 5 .0, LAST ) 

DO 24 1=1,12 

ITITLEt I ) “12 

I T I TLE( 1 ) “8HARDALAN - 

LABLE=4HG20 



CALL DRAW ( 500. X.G20. 2,0, LABLE, IT I TLE ,0*0, 0,0 ,0,0, 5, 5 .O.LAST ) 
DO 25 1=1,12 
25 I T I TLE ( 11*12 

I T I TLE ( 1 ) “8HARDALAN 
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LABLE=4HG30 

CALL DRAWI500»X»G30»2,0» LABLE • ITITLE»0»0»0»0»0»0»5»5»0»LAST ) 
DO 26 I « 1 * 12 

26 IT ITLE( I ) =12 

IT ITLE( 1 ) “8HARDALAN 
LABLE=4HG40 

CALL DRAW1500.X.G40, 3,0. LABLE, I TITLE *0,0. 0*0.0 ,0, 5. 5*0. LAST ) 
DO 27 1*1.12 

27 I T I TLE I I ) = 12 
ITITLEI 1 ) =8HARDALAN 
LABLE=4HC10 

CALL DRAW! 500.X.C10, 1.0. LABLE. I T I TLE *0 . 0 *0 * 0*0. 0.5. 5*0. LAST ) 
DO 28 1=1.12 

28 ITITLEC I ) *12 

I T I TLE ( 1 ) “8HARDALAN 
LABLE=4HC23 

CALL DRAW( 300. X.C23. 2.0. LABLE. I TITLE. 0.0. 0.0 ,0.0, 3. 5.0. LAST ) 
DO 29 1=1,12 

29 I T I TLE( I ) = 12 
ITITLEI 1 ) =8HARDALAN 
LABLE=4HC20 

CALL DRAW! 500 .X » C20 , 2 » 0 . LABLE ♦ I T I TLE .0 , 0 .0 ,0 ,0 .0, 5.5.0, LAST ) 
DO 30 1=1,12 

30 ITITLEC I ) >12 
ITITLEI 1)=8HARDALAN 
LABLE=4HC30 

CALL DRAW I 500,X,C30,2.0, LABLE , ITITLE, 0,0, 0.0.0, 0.5, 5,0, LAST ) 
DO 31 1=1,12 

31 I T I TLE I I ) = 12 
ITITLEI 1)=8HARDALAN 
LABLE=4HC40 

CALL DRAW I 5 00 .X.CAO. 3.0, LABLE, I T I T LE .0 . 0 .0 . 0 , 0 , 0 , 5 ♦ 5 .0 . LAST ) 
DO 32 1=1,12 

32 I T I TLE I I ) = 12 

I T I TLE 1 1 ) “8HARDALAN 
LABLE=4HL23 

CALL DRAW! 500.X.P23. 1.0 .LABLE, I T I T LE .0 , 0 .0 .0 , 0 ♦ 0 • 5 . 5 .0 , LAST ) 
DO 33 1=1.12 

33 ITITLEI I)=12 
ITITLEI 1)»8HARDALAN 
LABLE=4HL20 

CALL DRAWl 500. X, P20, 2 .0. LABLE. IT ITLE, 0,0,0, 0,0.0, 5, 5.0, LAST ) 
DO 34 1=1,12 

34 ITITLEI I ) =12 
ITITLEI 1)=8HARDALAN 
LABLE=4HL34 

CALL DRAW I 500. X.P34, 2,0. LABLE, IT ITLE.O, 0,0, 0,0, 0,5, 5,0 .LAST ) 
DO 35 1=1,12 

35 I T I TLE I I ) =12 
ITITLEI 1 ) =8HARDALAN 
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LABLE=4HL30 

CALL DRAW(500.X*P30.2»0.LABLE* I T I TLE *0 , 0 ,0 *0 *0 #0 * 5 * 5 *0*LAST ) 
DO 36 1=1*12 
36 I T I T L E ( 1 ) = 1 2 

1 T I TLE ( 1 ) =8 HARD ALAN 
LABLE=AHL40 

CALL DRAW(500*X.P40*3*0.LABLE* I T I T LE .0* 0*0 *0 .0 *0* 5 . 5 *0. LAST ) 

END 

END 
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Appendix III 

Computer Program for Example 3 . 1 



• • JOB0184F *ARDALAN 

PROGRAM NONLINER 

DIMENSION X( 500 ) »G1Q( 500) *G12 ( 500 ) *G20 ( 500 ) , CIO < 500 ) »C12 ( 500 ) .C20 ( 
1500) *ITITLE(12) 

A= 1 • 

GM12=1 • 

GM22*1 • 

G=l. 

0 = 1 . 

C2 = l* 

X( 1 )«0. 

DO 10 1=1*500 

G12(I)=G*EXPF(( GM22+GM12 ) *X ( I ) ) 

C12( I )=C2*(EXPF(GM22*X( I ) ) * ( EXPF ( GM12*X ( I ) )-l. ) ) 

C 1 0 ( I )=C1+C2*(EXPF(GM22*X( I ) ) * ( l.-EXPF ( GM12*X ( 1 ) ) )+ ( 1 .-EXPF ( GM12*X 
1(D) )**2) 

C20 ( I ) =-C12 ( I )+C2*EXPF(2.*GM22*X< I ) ) 

PRINT 15»X( I ) »G12 ( I ) *C10( I ) * C 1 2 ( I ) »C20( I ) 

15 FORMAT ( 1X*5F10.6) 

X( I+1)=X( I )+( A/500* ) 

10 CONTINUE 

DO 23 1=1.12 

23 I T I TLE ( I ) =12 

I T I TLE ( 1 ) =8 HARD ALAN 
LABLE=4HG12 

CALL DRAW ( 500 *X*G12 * 1.0.LABLE, IT ITLE.0.0.0. 0*0.0 *5*5 »0*LAST ) 

DO 24 1=1*12 

24 I T I TLE ( I ) * 12 

I T I TLE ( 1 ) =8HARDALAN 
LABLE=4HC10 

CALL DRAW ( 500.X* CIO* 2 *0*LABLE. I T I TLE *0* 0 *0 *0 *0 *0. 5* 5 *0* LAST ) 

DO 25 1=1.12 

25 ITITLEl I)=12 

I T I TLE ( 1 ) =8 HARD ALAN 
LABLE»4HC12 

CALL DRAW(500.X.C12.-2*0.LABLE* I T ITLE *0*0 *0 *0 .0 .0. 5* 5 .0* LAST ) 

DO 26 1=1*12 

26 I T I TLE( I ) *12 

I T I TLE ( 1 ) =8 HARD ALAN 
LABLE=4HC20 

CALL DRAW( 500 *X*C2Q* 3*0* LABLE. IT I TLE *0,0*0, 0.0 *0*5,5 *0,LAST ) 

END 

END 
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